Holographic analysis of diffraction structure factors 
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We combine the theory of inside-source/inside-detector x-ray fluorescence holography and Kossel 
lines/x ray standing waves in kinematic approximation to directly obtain the phases of the diffraction 
structure factors. The influence of Kossel lines and standing waves on holography is also discussed. 
We obtain partial phase determination from experimental data obtaining the sign of the real part 
of the structure factor for several reciprocal lattice vectors of a vanadium crystal. 
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The phase determination of diffracted beams is the 
central problem of x-ray crystallography. Several meth- 
ods exist to obtain this information, such as direct meth- 
ods pi, 3-beam diffraction M, anomalous diffraction, 
heavy atom or molecular replacement || and x-ray stand- 
ing waves or Kossel lines ||, |[ [| . Despite the many ad- 
vances in these methods, not all problems can be solved. 
Direct methods fail when the unit cell contains a large 
number of atoms. Anomalous diffraction and related 
methods are among the most successful methods for bi- 
ological crystallography, but they often require chemical 
modification of the molecules. Multiple beam diffraction, 
x-ray standing waves and Kossel lines have usually been 
applied only to high-quality crystals of relatively simple 
structures, or to the localization of dopants in high qual- 
ity crystals @. 

Kossel lines are formed when a source of short wave- 
length radiation (~ 1 A) is located on a crystallographic 
site: they result from the Bragg scattering of outgoing 
fluorescent x-rays from various sets of planes in the lat- 
tice. In the notation of holography, this is an "inside- 
source" experiment. The fine structure of these lines has 
been explained by the dynamical theory of x-ray diffrac- 
tion via the reciprocity theorem used in optics [4j. A 
proper analysis of the KL fine structure allows the de- 
termination of the phases of reflections j^, ||, 0, |l(], [ll| . 
In parallel to this work on KL, the x-ray standing wave 
(XSW) method has been developed §, | @, [jj. In 
this case, the source and the detector are interchanged 
as compared to the KL method: the atoms are subject 
to the changing wave-field in the crystal as the incident 
beam goes through a Bragg reflection, and fluorescent 
radiation proportional to the field at the atom is gener- 
ated. This constitutes the "inside-detector" configura- 
tion in holography [||. Besides the fine structure of the 
KL or XSW produced when the crystal orientation satis- 
fies a Bragg or Laue condition, tails are formed far from 
the Bragg angle. These coarse features, also formed by 
poorer-quality mosaic crystals provide information on the 
real and imaginary part of the structure factor jl0[ [LjJ . 

Unlike standard imaging methods, holography offers 



the possibility of extracting both intensity and phase 
information. X-ray fluorescence holography (XFH) is 
thus a very promising new method for obtaining a direct 
image in real space of the local environments of differ- 
ent atomic species in reasonably well-ordered crystals or 
molecular ensembles. Long-range translational order is 
not required, and indeed one of the most important re- 
sults obtained by XFH has been the imaging of the aver- 
age local environment of a quasicrystal [12] . Despite this, 
with the exception of quasicrystals, all systems measured 
so far have been well-ordered crystals lli 
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In this paper we develop a theoretical method for an- 
alyzing the inside-source and inside-detector holograms 
of a periodic object using the kinematic approximation. 
By analyzing the holographic reconstruction in reciprocal 
space, we discuss how to obtain directly the phase of the 
structure factors. We also show how the standard holo- 
graphic analysis is affected by diffraction, and discuss the 
solution to this problem. 

A hologram is formed whenever an unknown object 
wave E a bj is coherently added to a reference wave E re f: 



J(k) = \E Tei + E oh - 3 \ 
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where I re f and I bj are the intensities of the reference 
and object waves respectively. 

The hologram is extracted by subtracting the refer- 
ence beam intensity I re f and normalizing. The object 
term I i,j is usually assumed to be small and therefore 
neglected in order to permit a holographic analysis of 
the data in XFH. We will see later how this approxima- 
tion can affect our analysis. When an atom located at 
the origin emits (inside source) or detects (inside detec- 
tor) radiation, the resulting hologram (considering only 
the last term in eq. ([l]), i.e. the interference term) can 
be expressed as H: 



(k) = -2Rc^/(k,r) X r (k) 



(2) 
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where Xr (k) = — - — and /(k, r) is the scattering 



factor (including the Thompson scattering factor 17 ) of 
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the atom located at r. The summation is extended to all 
the atomic positions. 

Ideally, when the object term can be neglected 
and the scattering factors are isotropic, the holo- 
grams Xr (k) , Xr' (k) generated by two point scat- 
terers located at r and r', are orthogonal. i.e. 
(Xr (k) Xr' (k)) k —a (r) 8 (r — r'). a (r) is a normaliza- 
tion function and () k is the average on the measured 
k space. The two holograms are orthogonal because 
the relative phase oscillates when the infinite k-space 
is spanned. In normal situations, due in particular to 
the limited k-space sampling, 8 (r — r') becomes a point 
spread function peaked in r = r'. The typical reconstruc- 
tion algorithm is based on this orthogonality assumption. 
The holographic reconstruction is the projection of x (k) 
onto Xr' (k): 



C/(r') = (x(k)x:, (k)) t 



(3) 



What happens if we now perform the reconstruction tak- 
ing into account the long-range periodicity of the sys- 
tem? To answer this question, let us first examine the 
holographic reconstruction in the reciprocal space, by ap- 
plying a Fourier transform to eq. (|^). Since only the 
kernel Xr' (k) depends on r' we can bring the FT inside 
the average and write: 



G(h') = ^x(k) J d 3 re jh ' r X ;(k) 

= <x(k) X *,(k)) k , 
Xh (k) = / d 3 r^-e^ kr -^- h > r i = 



(|k-h| 2 -fe2) 



(4) 
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The formal analogy between (^) and (fij) suggests we 
consider Xh. (k) as the hologram generated by a point 
structure factor, located at position h of the reciprocal 
space, i.e. a sinusoidal charge density distribution with 
unitary scattering factor amplitude. Similarly to the dis- 
cussion above, we expect the holograms Xh' (k) , Xh (k) 
generated by two points h and h' in the reciprocal lat- 
tice to be orthogonal in k-space, i.e. (x h * (k) Xh (k)) — 
(3 (h) 8 (h — h') . In real situations the Dirac delta func- 
tion is replaced by a function peaked in h. 

We now examine the properties of the hologram x (k) 
in terms of the structure factors of the reciprocal lattice 
vectors. Following fTq] , we rewrite (||) in terms of the 
electron density p (r) : 



X(k) 



2Re Jd 3 xp(r) =^e 



i(kr— k-r) 



(G) 



where r e is the classical electron radius. For simplicity 
we have approximated the Thompson scattering factor 
to be a constant. If the system is periodic, the charge 
density distribution can be expressed in a Fourier series 
in terms of the reciprocal lattice vectors h and the rela- 
tive structure factors F^, namely — r e p (r) = F^e lh r , 



withF h = =pfp(r) 
v 



c d r, where V is the unit cell 



volume. The origin, which determines the phase shifts in 
the structure factors, is at the emitting atom in the unit 
cell. 

The hologram (|) becomes, in analogy with (0): 



x(k) = 2Re^F hXh (k) 



(7) 



The divergence of Xh (k) at the Bragg condition (h 2 — 2h- 
k = 0) in (||) is introduced as a consequence of neglecting 
extinction, mosaicity and the finite sample size. These 
effects can be approximated by writing k as a complex 
number [fill, i.e. k — k r +iki, and Xh (k) = Re{xh (k)} + 
Tm{xh (k)} then becomes: 



Re{xh (k)} 
Im{xh (k)} 
IXh(k)| 2 



|Xh (k)| 2 (h 2 -2h-k) 

|Xh (k)| 2 (-2^£v), 
l 



(8) 



As the wave- vector k changes across the pole of Xh (k) , 
there will be a symmetric and an antisymmetric con- 
tribution dependent in different ways on the real and 
the imaginary part of the structure factor. Writing 
Fh = Re {-Fh} + i Im {-Fk}, the hologram in (Q) becomes: 



X 



2^Re{xh}Re{F h }-Im{ Xh }Im{F h } . (9) 



Note that, at the Bragg condition, the object term can- 
not be neglected, it contributes with a second order diver- 
gence and therefore needs to be added to the interference 
term. The normalized object term can then be written 
as: 



Iob]_ 
Iref 



h 



l^h| 2 |xh| 2 



(10) 



Except when multiple Bragg conditions are satisfied si- 
multaneously, the second term can be neglected. We re- 
call that the object wave intensity is always smaller than 



that of the reference at a pole 
the interference term (| 



1 obj 
Iref 



< 1, while 



(2k r ki) 

at a Bragg condition, is 1 ™^' h . 
The interference term (|9j) can still be the dominant one, 
depending on the phase of the structure factor. 

Let us now consider the properties of the reconstructed 
image G, which can be viewed as the 'hologram of the 
reciprocal lattice', in relation to the structure factors. 
The real and imaginary parts of the reconstruction ker- 
ne l Xh' (k) (eq. rtq)) are respectively antisymmetric and 
symmetric across the pole as (h' 2 — 2h' • k) changes from 
positive to negative. The product between an antisym- 
metric and a symmetric function is antisymmetric, and 
the integral across the pole cancels out: 



(Rc{ Xh (k)}Re{x h (k)}) fc = a(h), 
(Re{xh(k)}Im{x h (k)}) fe » 0, 
(Im{xh(k)}Im{x h (k)}) fe = -6(h); 



(11) 
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where the normalization functions a (h) and b(h), are 
positive functions. The resulting holographic reconstruc- 
tion becomes: 

Re{G(h' = h)} » Re {F h } a (h) , 

Im{G(h' = h)} « ( 2 Im{F h }-g£)&(h). 

The real part of the reconstructed hologram, with h' 
equal to a reciprocal lattice vector h, is proportional to 
the real part of the structure factor. The characterization 
of Re{Fh} reduces the phase problem of the structure 
factor from the possible (0, 2tt) range to only two possible 
values of the phase, and in the centrosymmetric systems 
this correspond to the complete solution. This requires 
us to calculate the normalization function a (h), which is 
dependent on the measured k-space, polarization factors, 
extinction, sample boundaries, mosaic spread and angu- 
lar resolution. Such calculation goes beyond the scope of 
this work; however the knowledge of the sign of ReFh re- 
duces the range of the possible phases by half. The imag- 
inary part of G depends on the real and imaginary part of 
the structure factors. If the object term I D bj in ([l]) cannot 
be neglected, the Fourier transform of G, the "standard 
holographic reconstruction" , would have the wrong in- 
terpretation. It would be more some kind of Patterson 
function (the Fourier transform of -Fh| 2 ) rather than a 
hologram (Fourier transform of ReFh + ilmFh)- Only 
with the full knowledge of the normalization functions 
a (h) and b (h) can we separate and extract the values of 
the real and imaginary parts of the structure factors. The 
limitations on this straightforward analysis of the imag- 
inary part of the structure factor apply to the standard 
holographic reconstruction as well. It is not surprising 
that a diffraction measurement interpreted as a hologram 
in a recent article [ p0[ provides the correct structure, as 
the Patterson function peaks match the correct solution. 

Analyzing the reconstruction in the reciprocal space 
allows extracting directly information on the real part 
of the structure factors, simplifying the problem to the 
calculation of the normalization functions. By including 
the knowledge of the absolute value of the structure fac- 
tor, we can separate the contribution of the object inten- 
sity term to the imaginary part of the structure factor. 

We performed the experiments on a bending magnet 
beamline (BL 9.3.1) at the Advanced Light Source. A 
schematic drawing of the experimental setup is shown in 
Figure 1. Inside the chamber, the sample, a vanadium 
bcc single crystal with (111) orientation, is mounted on 
a standard two-axis goniometer. The V K a fluorescence 
radiation at 4.9 keV emitted by the sample was collected 
by a 4-channel high-speed solid state detector with single 
phot on p ulse analysis and a maximum of 4 MHz count 
rate |2l| . The measurement was performed by rotating 
the sample at high speed in a spiral motion (22|: the 
azimuth at 3600 degrees per second, and the polar at 
2 degrees per second and varying from (perpendicular 
to the surface) to 80 degrees. The detector was placed 
close to the sample to average the "inside source" holo- 




Detector 



FIG. 1: - Experimental setup. Monochromatic x-rays im- 
pinge on the sample, mounted on a two-axis goniometer that 
is rotated at high speed, and a solid-state detector collects 
the fluorescence radiation. Only the sample rotates, with the 
angle between the exciting x-rays and the detected fluorescent 
x-rays fixed. 




FIG. 2: - Measured inside-detector hologram of a vanadium 
(bcc) single crystal at 4.9 keV fluorescent energy and 6.3 keV 
incident energy, after symmetrization and rotation of the pat- 
tern to be in the (100) orientation. 



gram, resulting in an "inside-detector" measurement in 
holographic terminology The azimuth stepper mo- 
tor pulses were used to synchronize the data acquisition 
allowing us to collect a full pattern of 3.2 10 5 pixels in 
about 40 seconds. The measurement was repeated un- 
til the statistical noise and incident beam fluctuations 
have been reduced to a reasonable level; typically, sev- 
eral hundred separate patterns were thus summed in a 
final dataset. 

The measured hologram at 6.3 keV incident energy is 
shown in Figure 2 and the real part of the reconstructed 
'reciprocal' hologram Re{G(h')} as derived from Eqs. 
4 and 5 is shown in Figure 3. The network of white lines 
connects the known positions of the reciprocal lattice po- 
sitions in the V lattice. The image in reciprocal space 
shows its most intense yellow spots at the reciprocal lat- 
tice positions, which in turn correspond to the positive 
values of the real part of the structure factors. One can 
see that the height of these peaks is approximately con- 
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FIG. 3: Reconstructed structure factor obtained from the 
hologram in ^; axes scale is in 27r/A, the colorbar is in ar- 
bitrary units. Circles indicate the locations of the reciprocal 
vectors, with these being connected by light lines. At each 
circle, there is a positive peak, as expected for this system. 

stant. This is because the functions a (h) and b (h) are 
almost constant when h is not close to the wavenumber 
k, this can be true if the hologram is measured at a single 
energy in every direction, as was our case. The vanadium 
crystal is a special case in which the structure factors are 
real and positive. However the reconstructed image from 
a simulated Kossel line pattern of a PbSe single crystal 
which was chosen since it exhibit both positive and neg- 
ative signs in the structure factors, correctly shows these 
different signs. 

We have presented a method for the direct phase de- 



termination of the structure factors. This result has been 
obtained by analyzing the holographic reconstruction in 
reciprocal space and by combining the theory of insidc- 
source/inside-detector holography and Kossel lines/x-ray 
standing waves. This method can be applied to any crys- 
tal possessing an atom which can be excited to emit ra- 
diation. We have shown how holograms and standard 
holographic reconstruction can be distorted in periodic 
objects by x-ray diffraction, and discussed the possible 
solution to this problem. By separating the real and 
imaginary parts of the reconstructed image, and by calcu- 
lating the normalization function, we obtain the real part 
of the structure factor. By including the knowledge of the 
absolute value of the structure factor, one can separate 
the contribution of the object intensity term to imaginary 
part of the structure factor. 

We have demonstrated this method experimentally on 
a simple test case, by measuring the full inside-detector 
XSW pattern and obtaining the real part of the structure 
factors for a vanadium crystal. In order to obtain the full 
phase determination, the calculation of the normalization 
functions a (h) and b(h). However even the qualitative 
direct image obtained provides the sign of the real part 
of the structure factor, which for centro-symmetric sys- 
tems correspond to the complete solution. The informa- 
tion obtained by this technique can be used as input in 
a standard KL/XSW fitting analysis to obtain the full 
phase determination. 
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